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Recent measurements have shown oscillations in the upper critical field of simply connected meso- 
scopic superconductors. A quantitative theory of these effects is given here on the basis of a 
Ginzburg-Landau description. For small fields, the H — T phase boundary exhibits a cusp where 
the screening currents change sign for the first time thus defining a lower critical field H&. In the 
limit where many flux quanta are threading the sample, nucleation occurs at the boundary and the 
upper critical field becomes identical with the surface critical field H c z. 



I. INTRODUCTION 

As is well known, superconductivity only exists at suf- 
ficiently low temperatures T and small external mag- 
netic fields H . The resulting H — T boundary for the 
normal to superconducting transition is determined by 
the Ginzburg parameter k = A/£ . For type II supercon- 
ductors with k > l/y/2 one obtains a lower (H c \ ) and 
an upper ( H c % ) critical field which - for bulk samples 
- are universal functions of temperature [0J^]. However, 
it was realized long ago by Saint- James and de Gcnncs 
HQ that in the presence of a surface, these results are 
changed considerably. Regarding H ci , there is a sur- 
face barrier for the entrance of the first flux quantum. 
Thus the field up to which the sample stays in the Meiss- 
ner phase may be much larger than the thermodynamic 
H c i ||. In the case of the upper critical field, super- 
conductivity in a bounded sample persists even in the 
range H c2 < H < H c3 = 1.69H c2 provided the external 
field is parallel to the surface [Q. In this regime only a 
thin sheet at the sample boundary of the order of the 
zero field coherence length £ is superconducting. Quite 
generally, in samples whose size is of the order of the 
T = coherence length, one expects that the H—T phase 
boundary will strongly depend on the detailed form of the 
sample, reflecting the possible eigenmodes for the com- 
plex superconducting order parameter ijj(r) in the given 
geometry. Experimentally this was recently studied by 
Moshchalkov et al. , who investigated the temperature 
dependence of the upper critical field of small mesoscopic 
aluminium samples with typical sizes of lfi or less. The 
observed H — T phase boundary turned out to exhibit 
very peculiar size effects. Specifically it was found that 
Aharonov-Bohm like oscillations in the critical field were 
present even in a simply connected geometry, similar to 
the Little-Parks oscillations found long ago f| in thin 
walled superconducting cylinders. 

It is the purpose of this work to investigate size effects in 
the critical fields of small superconductors on the basis of 
a Ginzburg-Landau (GL) theory. We will find that such 
a description apparently remains valid down to system 



sizes of only a few coherence lengths. By a careful solu- 
tion of the boundary value problem for the linearized GL- 
equation near T c , we are able to quantitatively describe 
the observed structure in the upper critical field of a small 
disc. In the limit where many flux quanta are threading 
the sample, the upper critical field is in fact a surface crit- 
ical field, reproducing the standard H c ^ value of a semi- 
infinite geometry. Moreover we determine a generalized 
lower critical field for mesoscopic discs and rings. An 
interesting point is that the eigenvalue spectrum which 
determines the suppression of the critical temperature 
T C (H) as a function of magnetic field, is rather different 
from the case of electron levels in a quantum dot, because 
of the different boundary conditions. 

II. CRITICAL FIELDS OF SUPERCONDUCTING 
DISCS AND RINGS 

Let us consider a small disc with radius R and thick- 
ness d in an external magnetic field H = He z , which is 
perpendicular to the sample surface at z = id/2 . Near 
the normal to superconducting transition the change in 
the free energy with respect to the normal state can be 
expressed in terms of a GL-functional of the complex su- 
perconducting order parameter ip(r) 

m= ^ + /{£|(V-^A)?^(r)| 2 (1) 
v k 

+ ahMr)| 2 + !hMr)| 4 + f } d 3 r . 

Here B = V x A is the magnetic field in the sample with 
volume V, a = a' (T — T c )/T c and b are the standard 
GL-coefficients and [i the effective electron mass [pl^]. 
In principle there are also surface contributions to the 
free energy functional ([!]) which may be important for 
mesoscopic samples with a large surface to volume ratio. 
In our treatment below such contributions are neglected, 
which is justified only a posteriori. In the vicinity of the 
transition the order parameter and the screening currents 
are small. To lowest order we may therefore neglect the 



°Dedicated to Prof. W. Gotze on the occasion of his 60th 
birthday. 



1 



quartic term in F[ip] and replace the magnetic field by 
the external one. The most probable configuration of the 
order parameter which follows from the mean field equa- 
tion SFlipj/Stp* = is then determined by the eigenvalue 
problem 
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for a particle with charge 2e in an external magnetic field 
( e < ) . Assuming that the sample is embedded in an in- 
sulating medium, the relevant boundary condition is that 
of vanishing current normal to the sample surface dV . 
In covariant form the corresponding Neumann boundary 
condition is [fil 



zero temperature GL-coherence length, the maximum 
value of a - which always has m < - determines the 
magnetic field shift of the transition temperature by 



T c - T C (H) 
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with <l = 7r£ 2 (0)-ff and $ = hc/2\e\ the superconduct- 
ing flux quantum. In an infinite sample the ground 
state is the lowest Landau level with = Ticu c /2 , 
i.e. a°° = . The phase boundary is then given by 



T c - T C (H) 
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where n is a unit vector normal to the sample surface. 
In order to determine the H — T phase boundary, we 
must find the lowest eigenvalue E (H) associated with 
a nonzero order parameter ip(r) ^ 0. From Eq(H) the 
transition from the normal to the superconducting state 
is determined by 



Tc(H) 



E (H). 
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Here T c is the (mean field) transition temperature of the 
infinite system with zero field. Since Eq(H) > Eq(0) 
quite generally |7j], the transition temperature at finite 
field is always smaller or equal than at H = . In order 
to treat the case of discs or rings, it is convenient to intro- 
duce cylindrical coordinates (p, <f>, z) . In the appropri- 
ate gauge A = Hpe^/2 the solution of the Schrodinger 
equation (0) can then be written as 



t/> = K(p) e 



imcf> ^ik v z 
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Here m G is the angular momentum quantum num- 
ber and k v — vn/d with v £ ]N the discrete wavevector 
for motion in the z-direction. Since the lowest eigen- 
value has always v = , we will omit the z-dependence 
and the associated quantum number v in the follow- 
ing. Introducing a dimensionless variable ( = p 2 /2l 2 H 
with Ir = {hc/2\e\H) 1 / 2 the magnetic length for charge 
2e , the differential equation for H(C,) can be reduced to 
Kummer's confluent hypergeometric equation [pi 
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by the substitution TZ(Q — e 2 ^'2' w(() . The dimen- 
sionless parameter a is directly related to the eigenvalue 
E(H) by 
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with io c = \e\H/pc the standard cyclotron frequency. 
Using -Eq(H) = a and a' = h 2 /4p£ 2 {0) with £(0) the 



which is equivalent to the standard relation 
H = H c2 (T) = $ /2< 2 (T) |J|. For the finite system, 
the spectrum of eigenvalues follows from the Neumann 
boundary condition (|[) at the inner (Ri) and outer (R) 
radius of the ring. The general solution of (^) is a linear 
combination of Kummer functions B. In the case of a 
disc geometry only 



wi(C) = $(a, 



UHiJifaH + u) (10) 



is allowed since the second linear independent solution di- 
verges at the origin. Using standard recursion relations 
for 1F1 , it is straightforward to show that the bound- 
ary condition at p = R, which simply reads dR/dp = , 
since A • n = , leads to 

(|m| + l-a)$(a-l,|m| + l,Cfl) 
- $(a, |m| + 1, ( R ) + a $(a + 1, |m| + 1, ( R ) - (11) 

with (ji = R 2 /2lj I . For each given m equation jTl| ) 
determines a discrete series of eigenvalues a nm (H) , 
n € INq j which are decreasing with increasing n . They 
obey a nm < 1/2 and are continuous in ^ = <f>/$ |j] 
which is just the external flux $ = ttR 2 H through the 
area of the disc in units of the flux quantum. Analytical 
results for the spectrum can be obtained in the low field 
limit $ — * . In this limit it is straightforward to treat 
the general case of a ring with a = R4/R < 1 . Standard 
second order perturbation theory in the magnetic field 
then leads to a shift in the transition temperature which 
is given by 



T c - T C (H) = 1 
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The corrections to this result are of order $ 4 , since the 
ground state energy is even in $ . For a very thin ring 
with a — > 1~ this agrees with the low field limit of the 
Little-Parks result [p| 
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as expected. For general magnetic fields the phase 
boundary can only be obtained numerically. To this end 
we have directly solved the transcendental equation ( pd| ) 
which allows us to determine the spectrum without any 
discretization error. The energy levels are thus obtained 
with arbitrary accuracy, in contrast to previous work 
by Saint-James jlO| or by Nakamura and Thomas |Tl] ] 
who consider Dirichlet boundary conditions. The results 
are shown in Fig. 1, where the dimensionless eigenval- 
ues E nm for n — and m — 2, 1, 0, — 1, —10 are 
plotted as functions of &/$>o ■ 



their experiment, however, the presence of two gold con- 
tacts led to a boundary condition which is different from 
(^) over part of the sample boundary. While the os- 
cillations were still present, a detailed comparison with 
theory was difficult (for instance the first cusp was ob- 
served at $ pa 2.5$ compared to <I> = 1.924$ in the 
pure Neumann case). The more recent experiments of 
Moshchalkov et al. H, however, are in very good agree- 
ment with the theoretical predictions. This may be seen 
from a comparison with the measured deviation of the 
temperature shift AT C — T c — T C (H) from the average 
linear behaviour which is shown in Fig. 2. 




FIG. 1. Dimensionless energy eigenvalues of a disc as a 
function of the external magnetic flux. 



Evidently the lowest eigenvalue exhibits an oscillatory 
behaviour with cusps at values j = 1,2, ... where 

the magnetic quantum number of the lowest eigenstate 
jumps by one unit. The dimensionless distances between 
successive cusps 



<j>(i) _ 
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are given in table 1 with an accuracy corresponding to 
the last given digit. 



$W 1.923765 $ 


Ai = 1.923765 


$W = 3.392344 $ 


A 2 = 1.468579 


$(3) = 4.747920 $ 


A 3 = 1.355676 


$< 4 > : 6.045882 $ 


A 4 = 1.297962 


$ (5) 7.3068 $ 


A 5 = 1.2609 


$( 6 ) = 8.5423 $ 


A 6 = 1.2355 


$( y ) = 9.7584 $ 


A 7 = 1.2161 


$W = 10.9591 $ 


A 8 = 1.2007 


$< 9 > = 12.1477 $ 


A 9 = 1.1886 


$( 10 ) = 13.3255 $ 


Aio = 1.1778 



TABLE 1 

Experimentally the oscillatory behaviour of the ground 
state energy is directly reflected in the H — T phase 
boundary. For the case of a disc discussed here, this 
was actually first observed by Buisson et al. fll2fl. In 
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FIG. 2. Flux dependence of the oscillatory part of the tem- 
perature shift. The experimental data (squares) is taken from 
Ml, the solid line is the theoretical prediction based on equa- 
tion M). 



Here we have used the experimental value £(0) = 1/i , and 
a disc area which is only 2.7 % smaller than the area of 
the almost rectangular sample used in the experiment. It 
is important to note that the periods Aj decrease mono- 
tonically from Ai = 1.924 to Aoo = 1 (see table 1 and 
below) in contrast to an anomalous first period Ai w 1.8 
and constant successive ones A 2 « A3 ss A4 w 1.3 which 
were quoted by Moshchalkov et al. ||. 
The field at which the ground state changes from m = 
to m = — 1 allows us to extract a lower critical field 
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for a mesoscopic system with size R of order £ (0) . Here 
H c i is defined via the condition that for H < H c \ the 
sample tries to screen out the applied flux, whereas for 
H > H c \ the free energy is minimized by accepting one 
flux quantum. It is interesting to compare this with 
Fetter's theory of flux penetration in a superconduct- 
ing disc which is based on calculating the self en- 
ergy of a vortex. In the limit where the disc radius R 
is much smaller than the effective thin film penetration 
depth \2d — A 2 /d, it turns out, that it is energetically 
favourable for a vortex to enter if H > H c \ with fn3[ 
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Here r c sa £(0) is the core radius, which is always as- 
sumed to be much smaller than R. Obviously for samples 
whose size is of the order of the coherence length £ (0) , 
the expression Jig ) is no longer applicable. In this limit 
the approximation that the order parameter is constant 
beyond r c becomes invalid. As found above the lower 
critical field is then replaced by our result (|l5|), with a 
crossover at about R w 7£(0) . Here it is important that 
for R w £(0) , linearized GL-theory is sufficient to calcu- 
late Hci , because it is the sample boundary which lim- 
its the magnitude of the order parameter instead of the 
quartic term as usual. Finally consider a ring with inner 
radius R^ ^> r c . Then the lower critical field is simply 
determined by the condition that half a flux quantum is 
applied, i.e. 
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Indeed this follows from the quantization of the fiuxoid 
||, and is valid irrespective of the thickness of the ring. 
Comparing ( |l7| ) with the result jlif ) for a disc, we find 
that H c \ in the latter case is almost four times larger. 
Qualitatively this is due to the additional condensation 
energy in the center of the disc which is required for a 
vortex to enter. 

As a second point let us discuss the behaviour at <& 2> $0 
where many flux quanta have entered. In this limit the 
ground state has angular momentum |m| ^> 1 . The as- 
sociated eigenfunction is thus concentrated near the disc 
boundary. It is then obvious that our upper critical field 
for <I> 3> $0 is in fact a surface critical field. If this is 
correct, it should asymptotically approach the value ob- 
tained by Saint-James and de Gennes Q for a surface 
with a radius of curvature large compared to the coher- 
ence length. This can be verified by considering the spe- 
cial values $ m , m = 1, 2, ... in Fig. 1, where the tangent 
to Eo\ m \($) goes through the origin (i.e. we are consid- 
ering successive approximations to the envelope). These 
values are given in table 2 together with the correspond- 
ing values of a max . 
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$/$ 




T C -T C {H) 1 $ 

T / 

± r SPO 


1 


2.44 


0.28761 


0.849 


2 


3.92 


0.26561 


0.937 


3 


5.28 


0.25514 


0.979 


4 


6.56 


0.24872 


1.005 


5 


7.82 


0.24426 


1.023 


6 


9.09 


0.24091 


1.036 


7 


10.28 


0.23839 


1.046 


8 


11.46 


0.23638 


1.054 


9 


12.68 


0.23449 


1.062 


10 


13.86 


0.23300 


1.067 


100 


110.23 


0.21395 


1.144 


200 


214.72 


0.21130 


1.154 


1000 


1033.76 


0.20778 


1.169 


10000 


10109.33 


0.20583 


1.177 



TABLE 2 

It is obvious that a ma x converges to a limiting value 
0.2058... . Using (g) the associated transition tempera- 
ture is then given by 



T c - TJH) $ 



(18) 



As expected this is completely equivalent to the well 
known result H = H c3 (T) = 1.695ff c2 (T) for the surface 
critical field Q . The coefficient 1 — 2a™ aa , < 1 is in fact 
just the ratio between the ground state energy in the 
disc and the energy h,u> c /2 of the lowest Landau level in 
an infinite sample. Edge states centered near the sam- 
ple boundary have thus a lower energy than bulk levels. 
Note that this behaviour is just the opposite of the case 
with Dirichlet boundary conditions (relevant e.g. for edge 
states in the Quantum Hall Effect) where ed ge s tates are 
above the corresponding bulk Landau levels jjjj . Finally 
let us discuss the behaviour of the periods of the 
ground state oscillations for large flux $ ^> $0 • Due to 
the factor pl m l/ 2 the order parameter for increasing mag- 
netic quantum number \m\ 3> 1 is more and more con- 
centrated near the sample boundary, but is practically 
zero in the interior of the sample. The simply connected 
disc thus effectively behaves like a ring with a normal 
core of size R — c±Ih , where c\ is a constant Jll|. The 
periodicity observed in Eq (H) is then simply determined 
by the condition that one additional flux quantum enters 
the area of the normal core, i.e. 



Aj($ > $ ) = 1 



2 Cl ^ 
R 



1. 



(19) 



In fact this field dependance was observed in the exper- 
iments by Buisson et al. already for $/$o > 5 (l2). In 
the asymptotic limit, which is however only reached for 
<&/<l>o > 10 3 (see table 2) the coefficient c\ can be ob- 
tained analytically as 2ci = \/0.59 ~ 0.76 p0[|. 
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III. DISCUSSION 

Using linearized GL-theory we have calculated the nu- 
cleation field of a small superconducting disc with a ra- 
dius which is of the order of the coherence length £ (0) . 
The good agreement with the experimentally observed 
H — T phase boundary suggests that the macroscopic 
GL-description remains valid in this regime which is not 
obvious a priori. A surprising feature of our results is 
that Aharonov-Bohm like oscillations are present even in 
a simply connected sample. The physical origin of this ef- 
fect is that already the entrance of a single flux quantum 
effectively makes the sample a multiply connected one. In 
the limit $ ^> $o the disc behaves like a thin walled ring, 
leading to oscillations in T C (H) which are completely 
equivalent to the well known Little-Parks experiment. 
It is interesting to note that these effects depend cru- 
cially on the Neumann boundary conditions. In fact the 
equivalent eigenvalue spectrum with Dirichlct boundary 
conditions, which was studied by Nakamura and Thomas 
Pt , does not exhibit any oscillations in the ground state 
energy E ($) . It is an interesting future problem to in- 
vestigate similar effects in the fluctuation diamagnetism 
or extend the calculations above to more complicated 
geometries. This would allow to study eigenvalue spec- 
tra for systems with classical chaotic dynamics [ p^|Jl6| 
without the complications due to electron-electron in- 
teractions which are unavoidable in non-superconducting 
mesoscopic systems. 
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